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ABSTRACT 

Context. Probing the structures of stellar winds is of prime importance for the understanding of massive stars. Based on their optical 
spectral morphology and variability, the stars of the Oef class have been suggested to feature large-scale structures in their wind. 
Aims. High-resolution X-ray spectroscopy and time-series of X-ray observations of presumably-single 0-type stars can help us un¬ 
derstand the physics of their stellar winds. 

Methods. We have collected XMM-Newton observations and coordinated optical spectroscopy of the 06 lef star A Cep to study its 
X-ray and optical variability and to analyse its high-resolution X-ray spectrum. We investigate the line profile variability of the He ii 
A 4686 and Ha emission lines in our time series of optical spectra, including a search for periodicities. We further discuss the variabil¬ 
ity of the broadband X-ray flux and analyse the high-resolution spectrum of A Cep using line-by-line fits as well as a code designed 
to fit the full high-resolution X-ray spectrum consistently. 

Results. During our observing campaign, the Hen 44686 line varies on a timescale of ~ 18hours. On the contrary, the Ha line 
profile displays a modulation on a timescale of 4.1 days which is likely the rotation period of the star. The X-ray flux varies on time- 
scales of days and could in fact be modulated by the same 4.1 days period as Ha, although both variations are shifted in phase. The 
high-resolution X-ray spectrum reveals broad and skewed emission lines as expected for the X-ray emission from a distribution of 
wind-embedded shocks. Most of the X-ray emission arises within less than 2R, above the photosphere. 

Conclusions. The properties of the X-ray emission of A Cep generally agree with the expectations of the wind-embedded shock 
model. There is mounting evidence for the existence of large-scale structures that modulate the Ha line and about 10% of the X-ray 
emission of A Cep. 
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1. Introduction 

High-resolution X-ray spectroscopy with the RGS and the 
HETG instruments onboard of XMM-Newton and Chandra, re¬ 
spectively, has revolutionized our view of the X-ray emission of 
massive stars. One of the first O-type stars that was observed 
in X-rays at high spectral resolution was the 04 lef supergiant 
^ Pup (Cassinelli et al. l2001l Kahn et al. l2001l) . often considered 
as the prototype O-supergiant. The spectrum revealed broad, 
skewed, and blue-shifted lines, as expected for an X-ray emitting 
plasma distributed throughout the stellar wind and moving along 
with this wind (MacEarlane et al. ll991l Owocki & Cohen l20011) . 
However, a detailed analysis of the shape of the spectral lines in¬ 
dicated that there is less absorption in the wind than expected 
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given the mass loss rate of the star inferred under the assump¬ 
tion of a homogeneous wind. Moreover, observations of a small 
sample of other, relatively nearby and hence X-ray bright, single 
O-type stars with XMM-Newton and Chandra revealed lines that 
were less broad, less skewed and less blue-shifted than expected 
(see Giidel & Naze 120091 and references therein). These findings 
were interpreted as evidence for the winds of O-type stars be¬ 
ing clumpy (or even porous) and the mass-loss rates being lower 
than previously thought (Eeldmeier et al. l20031 Owocki & Cohen 
120061 Oskinova et al. l2^ Cohen et al. l20T0llm^ . 

( Pup displays X-ray line morphologies as predicted for a 
wind-embedded plasma, but has several properties that make it 
stand out from the sample of presumably single O-stars observed 
at high-resolution in X-rays. For instance, it has the largest mass- 
loss rate and one of the earliest spectral types among the stars of 
this sample. Moreover, its optical spectrum displays peculiarities 
that are partly reflected by its Oef spectral classification. 

In this context, we decided to observe another Oef star, as 
similar as possible to ( Pup. The best candidate for this purpose 
was the 06 lef supergiant 4 Cep (= HD 210839). The peculiar¬ 
ities of the optical spectrum of 4 Cep and its similarities with 
( Pup were already noted by Wilhelmina Fleming (Pickering & 
Fleming |T905]|. In the early 1970’s, it was found that both 4 Cep 
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and { Pup belong to a handful of stars that display a double- 
peaked He II /14686 emission line in their spectrum. These ob¬ 
jects were classified as so-called Onfp stars by Walborn (I1973I I 
or Oef stars by Conti & Leep (I1974I I. ^ Pup (V = 2.25) and A Cep 
(V = 5.05) are by far the brightest members of this class (both 
in optical and X-rays). In addition, both ( Pup and A Cep are 
most-probably single runaway stars (Gies I1987t and are rather 
fast rotators (Conti & Ebbets \\911\ Penny [199^ Howarth et al. 
119971 Simon-Diaz & Herrero 12007112014I I. The peculiar mor¬ 
phology of the He ii A 4686 line of Oef stars is likely related to 
this fast rotation. Indeed, Hillier et al. (I2012I I demonstrated that 
the double-peaked profile of this line stems from a complex in¬ 
terplay of emission and absorption of a rotating wincQ. Hillier et 
al. (I2012I I further showed that the box-shaped morphology of the 
N III AA 4634-4642 blend is also due to the effect of rotation, but 
this time on a purely photospheric emission line. Harries (120001 1. 
Harries et al. (I2002I I and Vink et al. (I2009I I reported complex lin¬ 
ear polarization effects across the Ha line profiles of (, Pup and 
A Cep. These authors attributed these effects to the impact of the 
rapid rotation of the stars on their winds. 

The idea of a corotating wind coupled with presumably ro¬ 
tational modulations of the optical emission lines of ^ Pup led 
Moffat & Michaud (I1981I I to the suggestion that the winds of 
Oef stars might be shaped by a dipolar stellar magnetic field. 
However, David-Uraz et al. (I2014I I recently revisited the mag¬ 
netic properties of a sample of OB stars, including ^ Pup and 
A Cep. Their goal was to establish whether or not the discrete 
absorption components (DACs) seen in the UV spectra of their 
sample stars are associated with large-scale dipolar fields. No 
significant longitudinal magnetic field was observed in any of 
the stars; for A Cep an upper limit on the dipolar field strength 
of 136 G was inferred, corresponding to a wind confinement pa- 

rameter 77 , = * <0.15 which should have no impact on the 

global geometry of the wind of A Cep nor on its X-ray emission. 

In this paper, we present the very first analysis of the X-ray 
spectrum and lightcurve of A Cep, as well as the results of a 
simultaneous monitoring of the star in optical spectroscopy. Our 
observations are presented in Sect.|2 The variability of A Cep, in 
the optical and the X-ray domain, is addressed in Sects. [3]and|4] 
respectively. Sect.|5]presents our analysis of the X-ray spectrum 
of A Cep. Finally, Sect.| 6 ]discusses our results and presents our 
conclusions. 


2. Observations 

We collected four observations of A Cep with XMM-Newton 
(Jansen et al. 1200H . Each of these observations had a nominal 
duration of 70 - 95 ks (see Table[TJ. Observations I and III were 
affected by a rising background level towards the end of the ob¬ 
servation (typically during the last 10 - 15% of the exposure). 
The corresponding time intervals were discarded. 

The EPIC instruments (Turner et al. 120011 Striider et al. 
12001b were operated in full frame mode and with the thick filter 
to reject optical and UV photons. The thick filters efficiently sup¬ 
press the optical contamination for sources with V > 4 (MOS 1 
& 2) and V > 1 (pn). Given the optical magnitude of A Cep 
(V = 5.05) we thus expect no optical loading of the EPIC 
CCDs. The reflection grating spectrometer (RGS, den Herder et 
al. 12001b was operated in normal spectroscopy mode. 


* Hillier et al. l l2012b assume a rigidly rotating wind out to a distance 
where the wind velocity reaches 20 km s“', and conservation of angular 
momentum beyond that point. 


Table 1. Journal of the XMM-Newton observations 


Obs 

Date 

JD-2450000 

Exposure time 
(ks) 

I 

6456.677 

82.5 

II 

6508.319 

75.8 

III 

6510.473 

94.9 

IV 

6514.317 

70.7 


Notes. The date of the observation is given at mid-exposure. The expo¬ 
sure time indicated in the third column corresponds to the RGS 1 instru¬ 
ment. 


The raw data were processed with the Scientific Analysis 
System (SAS) software version 13.0. The EPIC spectra of the 
source were extracted over a circular region centered on the 
2 MASS coordinates of A Cep and adopting an extraction radius 
of 50arcsec. The background spectrum was evaluated from a 
nearby, source-free circular region located on the same detec¬ 
tor chip as the source itself. The RGS 1 and RGS2 first and sec¬ 
ond order spectra were extracted for each observation individu¬ 
ally, and were then combined for the subsequent analysis. Since 
A Cep is a rather isolated X-ray source, there is no risk of con¬ 
tamination of its RGS spectrum by nearby bright X-ray emitters. 
Lightcurves of the source and the background were extracted 
for each EPIC instrument and converted into equivalent on-axis, 
full-PSF count rates using the epiclccorr command under SAS. 
We built lightcurves over the full energy band and for three dif¬ 
ferent time bins (100 s, 500 s and 1 ks), as well as over soft (0.2 
- l.OkeV) and hard (1.0 - lO.OkeV) bands with 1 ks time bins. 

In support of the X-ray observations, we collected optical 
spectroscopy. For this purpose, we used the AUITELIE spec¬ 
trograph (Gillet et al. 11994b at the 1.52 m telescope of the 
Observatoire de Haute Provence (OHP, France) and the refur¬ 
bished HERDS spectrograph at the TIGRE telescope (Schmitt 
et al. 12014b at La Luz Observatory (Guanajuato, Mexico). The 
AURELIE spectra have a resolving power of 7000 and cover 
the region between 4450 and 4880 A. The exposure times were 
10 min per observation and the mean S/N ratio was 340. Between 
two and four consecutive exposures were taken per night. The 
HERDS echelle spectra have a resolving power of 20 000 and 
cover the full optical range, though with a small gap near 5800 A. 
The integration times were 5 min, resulting in a mean S/N of 
160 per pixel. The spectra were usually taken by series of three 
consecutive integrations. Whenever possible, several series were 
taken during the night (beginning, middle and end of the night). 

A first set of optical spectra was taken in June 2013 (16 
AURELIE spectra and 3 HERDS spectra) around the time of X- 
ray observation I. A more intensive monitoring was performed 
in August 2013 covering the X-ray observations II to IV (35 
HERDS spectra). Figure[T] illustrates the EPIC-pn lightcurve of 
A Cep along with the sampling of the optical spectra. 

The HERDS data were reduced with the corresponding 
reduction pipeline (Mittag et al. 120101 Schmitt et al. 12014b 
whilst the AURELIE spectra were processed with the MIDAS 
software provided by ESD. The spectra were normalized 
self-consistently using a series of continuum windows. 

Unless otherwise stated, the quoted error bars in our analysis 
of the optical and X-ray data refer to 1 cr errors. 
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Fig. 1. The EPIC-pn lightcurve of A Cep is shown along with 
tickmarks that indicate the times of the AURELIE and HEROS 
spectra. The top and bottom panels correspond to the June and 
August 2013 campaigns, respectively. 


3. Optical Variability 

Spectroscopic variability in the optical and UV domain is 
widespread among O supergiants. This variability can arise from 
small- and large-scale structures in the stellar winds, non-radial 
pulsations, or magnetic helds. In the case of A Cep, the vari¬ 
ability of the He ii A 4686 and Ha emission lines was hrst noted 
by Conti & Erost (119741 1. The variability of the star’s emission 
lines has been subsequently studied by a number of authors 
(Lacv [T977] keep & Conti [T^ Ebbets [T982l Grady et al. [19831 
McCandliss 119881 Henrichs et al. 119881 Henrichs ll9911 Kaper 
et al. 119971119991 Markova et al. 120051) . Although some quasi- 
cyclic behaviour with a typical timescale near 2 days, interpreted 
as half the rotation period, was found at some epochs (Henrichs 
& Sudnik 12014t . the variations do not seem to be strictly 
(mono)periodic, but are reminiscent of the epoch-dependent 
variability found in other Oef stars (Rauw et al. 120031 De Becker 
& Rauw [2004lFl Henrichs & Sudnik (120141 1 accordingly suggest 
that the variations might arise from the presence of several tran¬ 
sient short-lived magnetic loops corotating with the star. The lo¬ 
calized magnetic helds responsible for these loops could take 
their origin in sub-surface convective layers of the star. 

Absorption lines were also found to be variable (e.g. 
Eullerton et al. 119961) . and de Jong et al. (119991 see also 
Kholtygin et al. 1201 111 claimed the detection of low-order non- 
radial pulsations (NRPs) with periods of 12.3 (/ = 3) and 6.6 hrs 
(/ = 5). However, a recent multi-epoch analysis of the absorption 


^ Quite recently, Howarth & Stevens ( 120141 reported however on a 
period of 1.78 days in spacebome photometry of { Pup. This period 
was found to be stable over the three year interval covered by the data 
and was attributed by Howarth & Stevens ( 120141 to low-i pulsations. 


line prohle variability revealed a more complex situation (Uuh- 
Sonda et al. 120141120151 . The frequency content of the power 
spectrum of these spectral variations considerably changes from 
one epoch to the other and no stable frequency was found that 
could unambiguously be attributed to pulsations. As becomes 
clear from Eig.[T] the sampling of our new optical data is not 
suited to study such rapid variations. 

In the study of our optical time series, we accordingly focus 
on the longer-term variability exhibited by two emission lines 
(He II A 4686 and Ha) which are at least partially formed in the 
wind. The Ha line is heavily affected by telluric absorption fea¬ 
tures that need to be corrected to avoid biasing our variability 
analysis. As a first step, we thus used the telluric tool within 
IRAE along with the list of telluric lines of Hinkle et al. (120001 ) 
to remove these absorptions in the Ha region. 

The time-series of continuum-normalized and telluric- 
corrected spectra were analysed by means of dedicated tools. 
The spectroscopic variability was first quantified via the tempo¬ 
ral variance spectrum (TVS, Eullerton et al. I1996t . In the blue 
spectral range, the strongest variability affects the He ii A 4686 
emission line. Beside the Hei 44471, Heii 44542, and Niii 
44 4509-4522 absorptions, and the N iii 44 4634-4642 emission 
lines discussed by Uuh-Sonda et al. (120141 ). strong variability is 
also found in the Hj3 absorption line. In the red spectra, the Ha 
emission displays by far the strongest line-profile variability. 

To characterize the timescales of the line profile variability, 
we have then applied a 2D-Eourier analysis (Rauw et al. l2008t . 
In this method, the Eourier periodogram is computed at each 
wavelength bin by means of the algorithm of Heck et al. (119851 ) 
modified by Gosset et al. (120011 ). This latter method explicitly 
accounts for the irregular sampling of the time-series. 

Eor the August 2013 dataset, we hnd that the periodogram 
of the He II 44686 line displays a strong peak at vi = (1.315 + 
0.014)day“' (or P = 18.3 + 0.2hrs). We have prewhitened 
the data for the sinusoidal signal corresponding to this fre¬ 
quency. Eigure|2] illustrates the mean periodogram before and 
after prewhitening. Erom this figure, one clearly sees that vi ac¬ 
counts for most of the power in the periodogram. The changes 
of the amplitude and phase constant of the V] modulation as a 
function of radial velocity are shown in Eig.|3] The errors were 
estimated using Monte-Carlo simulations (see Rauw et al. 120081 ). 
The modulation extends over the full width of the line, although 
its amplitude is largest in the central reversal between the two 
emission peaks. 

The other lines in the blue domain that display significant 
variability, the absorption lines He 1 4 4471, He ii 4 4542, and Hj8 
have distinct periodograms that are dominated by two frequen¬ 
cies, 0.46 and 0.57 day”', which are aliases of each other. This 
points towards a different mechanism for the variations of photo- 
spheric and wind lines. We stress here that the detected frequen¬ 
cies should not be considered as reflecting strict periodicities. 
Indeed, 4 Cep and other Oef stars are known to display vari¬ 
ability with timescales that can change quite significantly from 
one epoch to the other (see above). Moreover, some part of the 
variability of the emission lines is likely due to statistical fluctu¬ 
ations of the number of small-scale clumps in the wind (see e.g. 
Eversberg et al. 119981) and is thus not expected to follow a peri¬ 
odic behaviour. Our current results should therefore be regarded 
as indications of the timescale of variations at the epoch of the 
X-ray observations. 

The periodogram of the August 2013 Ha data (see Eig.lH 
reveals no evidence for modulation at the vi frequency. Instead, 
most of the power in the periodogram is due to two frequencies 
V 2 = (0.245+0.014)day-' andvj = (0.495±0.014)day-'. The V 3 
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Fig. 2. Periodogram of the He ii A 4686 emission line as observed 
in August 2013. In the mean periodogram over the wavelength 
interval 4678 - 4692 A, as computed from the time series of 
HEROS spectra (top panel), the highest peak is found at the 
frequency vi = 1.315 day“^. The residual periodogram, after 
prewhitening the data for the signal corresponding to the vi fre¬ 
quency, is shown in the middle panel. The bottom panel provides 
the power spectral window of our time series. 
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Fig. 3. Variations of the Heii 44686 emission line with the vi 
frequency in August 2013. From top to bottom, the panels illus¬ 
trate the mean line profile, the amplitude and the phase constant 
of the modulation as a function of wavelength. 


frequency is probably the first overtone of V 2 , but there does not 
seem to exist any aliasing or other logical link between vj and 
any of V 2 or V 3 . It is remarkable here that V 2 is fully compatible 
with the suspected rotation frequency of A Cep near 0.25 day ' 
(see Henrichs & Sudnik l20141 l. although, given the unstable be¬ 
haviour of the variability of A Cep, this could also be a coinci¬ 
dence. An interesting point is that the V 3 frequency is quite close 
to the 0.46 day“' frequency seen in the power spectra of the ab¬ 
sorption lines in the blue part of the spectrum (see above). This 


Fig. 4. Same as Fig.|3 but for the Ho- emission line as observed 
in August 2013. The mean periodogram over the wavelength 
interval 6545 - 6580 A (top panel) reveals its highest peak at 
V 2 = 0.245 day '. The middle panel displays the residual pe¬ 
riodogram once the data have been prewhitened for the signal 
corresponding to the V 2 and V 3 frequencies. 
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Fig. 5. Variations of the Ha emission line with the V 2 and V 3 fre¬ 
quencies in August 2013. From top to bottom, the panels illus¬ 
trate the mean line profile, the amplitudes (black and red sym¬ 
bols stand for the variations at frequency V2 and V3, respectively), 
and the phase constants of the modulation as a function of wave¬ 
length. 


suggests that part of the variability of Ha takes its origin in the 
same region as the absorption lines. 

Assuming that V 3 is indeed the n - 2 harmonic frequency 
of V 2 , we have used the information on phase constant and am¬ 
plitude as a function of wavelength (Fig.|5]l to reconstruct the 
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Fig. 6. Variations of the EW of the Ha line. Phase 0.0 corre¬ 
sponds to the time of the first optical spectrum of the August 
2013 campaign (HJD 2456507.7607). The solid curve shows 
the expected variations as reconstructed from the amplitudes and 
phases of the V 2 and 2xv2 frequencies (see Fig.|5j, whilst the data 
points with error bars correspond to the actual EWs measured in 
August 2013. 


variations of the Ha line within the cycle corresponding to V 2 . 
We have then computed the equivalent width (EW) of the recon¬ 
structed profiles as a function of phase. The comparison with the 
observed EWs is shown in Eig.|6l Eor the observations, the errors 
were evaluated using the formalism of Vollmann & Eversberg 
(I2006t . Although the dispersion is quite large, the dominant 
trends of the data are indeed predicted by the results of the re¬ 
construction. 

Due to poor weather conditions at the TIGRE site, the June 
2013 dataset is less densely sampled. The OHP data do not cover 
the Ha line and we thus restrict the analysis of the June 2013 
data to the He ii A 4686 line. Because of the poorer sampling, the 
Fourier analysis is subject to a much more severe aliasing. The 
highest peak at (1.405 + 0.02) day“^ accounts again for most of 
the power. Although it is close to vi, the severe aliasing leads to 
a strong ambiguity about the actual frequency of the signal. 

Our hnding of a very different behaviour of the He ii A 4686 
and Ha lines of A Cep is surprising at hrst sight, as it would 
suggest that these emission lines arise from different parts of 
the wind. Under local thermodynamical equilibrium (LTE) con¬ 
ditions, one would expect the Hen 44686 to form closer to 
the photosphere than Ha given the very different excitation po¬ 
tentials of the upper levels of the He ii and the Ha transitions 
(50.8 versus 12.1 eV). However, the atmospheres of massive 
stars strongly deviate from LTE. For instance, the populations 
of the energy levels involved in the He ii A 4686 transition are 
affected by pumping due to the He ii Lyman lines below 300 A. 
Figure|2] illustrates the contribution functions ^ of both lines as 
a function of height above the photosphere as derived from the 
best-fit non-LTE CMFGEN model of Bouret et al. (120121 1. The 
contribution function f is dehned by 

^OO 

EW - I ^dlogr 
Jr, 



r/R. - 1.0 


Fig. 7. Normalized contribution functions showing the location 
of the line formation regions of the He ii A 4686 (short-dashed 
line) and Ha (long-dashed line) transitions according to the best- 
fit CMFGEN model of Bouret et al. (120121 1. The solid line pro¬ 
vides the corresponding velocity law. 


(see Hillier [T9891 l. For the purpose of the discussion here, it is 
sufficient to consider the relative spatial variations of In Fig. 
|7] we have thus normalized the contribution functions of both 
lines by dividing by their maximum value. In this spherically 
symmetric model, both lines actually form over similar regions, 
although the photospheric contribution to the He ii A 4686 line 
arises deeper in the photosphere and this line’s contribution func¬ 
tion drops slightly faster in the wind than the one of Ha. Whether 
or not this slight difference is sufficient to explain the differences 
in the observed behaviour is unclear. In this context, Martins et 
al. (I20I5I) recently studied the radial dependence of variability 
in seven 09 - BO.5 stars. They found that lines that form in the 
photosphere can show very different variability patterns despite 
their line formation regions being close to each other. Our study 
indicates that this conclusion also applies to the case of the much 
earlier 06 lef star A Cep. 

In addition to the radial stratihcation discussed above, one 
might wonder whether the fact that A Cep is a fast rotator could 
lead to a latitudinal separation of the line formation regions. 
Indeed, in fast rotators such as HD 93 521, gravity darkening 
is sufficient to produce a latitudinal separation of the formation 
regions of low-excitation absorption lines, which form near the 
cooler equator, and high-excitation lines, which are preferen¬ 
tially formed near the hotter polar regions (Rauw et al. l2012l and 
references therein). Adopting the Roche model, a rotational pe¬ 
riod of 4.1 days and the stellar parameters (R,, M») from Bouret 
et al. (120121 1. we estimate that for A Cep the effective temper¬ 
atures near the pole and near the equator would differ by 3%, 
i.e. by about 1200 K. This effect might also contribute to the ob¬ 
served difference in variability behaviour of the He ii A 4686 and 
Ha lines, although a decrease in temperature of about 5000 K 
would be required to reach the temperature regime where the 
He II line is most sensitive to temperature. 
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4. X-ray lightcurves 

Since X-ray emission from single non-magnetic O-type stars 
is commonly considered to arise from a distribution of shocks 
inside the stellar wind, one expects this X-ray emission to be 
variable on timescales of hours on which these shocks evolve. 
However, a comprehensive search for such short-term variabil¬ 
ity in the case of f Pup provided no signihcant detection (Naze 
et al. I2013II . This result implies a lower limit on the number of 
X-ray emitting parcels in the wind of f Pup of 10^ (Naze et al. 
12013b . Still, Naze et al. (12013b found significant trends in the X- 
ray emission level on longer timescales of at least several days. 
These variations could be related to large-scale structures ex¬ 
tending over a significant fraction of the wind and corotating 
with the star. In this respect, it is interesting to look for corre¬ 
lations between the optical and X-ray variability. Such a corre¬ 
lation between the Ha emission strength and the ROSAT X-ray 
count rate in a narrow energy band had previously been reported 
by Berghofer et al. (11996b for ^ Pup. The data from both wave¬ 
length ranges revealed a period of 16.7 hrs. However, subsequent 
studies with ASCA (Oskinova et al. 12001b and XMM-Newton 
(Kahn et al. 120011 Naze et al. 12013b did not conhrm the pres¬ 
ence of this periodicity, suggesting that it was probably a tran¬ 
sient feature. Another important result in this context was the 
recent discovery by Oskinova et al. (12014b of X-ray pulsations 
in the fi Cep-like pulsator CMa. The X-ray pulsations occur 
on the same period as the optical pulsations suggesting a direct 
connection between the two phenomena. Since both pulsations 
and large-scale wind structures have been reported for A Cep, we 
decided to take advantage of our data to study the variations of 
its X-ray flux. 

We started by investigating the intra-pointing variability (see 
Fig.|8]l. Following the approach of Naze et al. (12013b . we per¬ 
formed tests of the EPIC-MOS and pn background- and 
barycentric-corrected lightcurves of each observation. We did so 
for time-bins of 100s, 500s and Iks, keeping only bins with 
more than 80% fractional exposure time, and for several hy¬ 
potheses (constancy, linear trend, quadratic trend). No signifi¬ 
cant intra-pointing variability is detected with these tests for ob¬ 
servations I, II and IV. Significant, i.e. with significance level 
(SL) < 1%, intra-pointing variability is found for the pn data 
and bins of 100 s in observation III. The MOS data of this ob¬ 
servation have larger SL, hence show only moderate variabil¬ 
ity. However, for the 100 s bins, the relative dispersions of the 
count rates about the mean are only marginally larger than what 
is expected from Poissonian photon noise. Some datasets appar¬ 
ently exhibit oscillations of the count rate. This is the case for 
instance of the pn data of obs. II (see Fig.|8]l where a Fourier 
analysis yields a frequency of 4.48d '. However, this signal is 
not present in the MOS data and does not exist in the pn data of 
the other observations. 


Table 2. EPIC count rates of A Cep in the 0.2 - lOkeV energy 
band averaged over the different observations. 


Obs. 


Count rates (cts s ‘ 



MOSl 

MOS2 

pn 

I 

0.117 ±0.001 

0.119 ±0.001 

0.430 ± 0.003 

II 

0.124 ±0.001 

0.124 ±0.001 

0.455 ± 0.003 

III 

0.125 ±0.001 

0.126 ±0.001 

0.459 ± 0.003 

IV 

0.126 ±0.002 

0.129 ±0.002 

0.466 ± 0.003 


Notes. The Icr uncertainties of the mean are quoted. 






Time (ksec) 


Fig. 8. Background-corrected EPIC-pn X-ray lightcurves of 
A Cep with a time-bin of 1 ks. The abscissae show the time 
elapsed since the start of each observation. 


Since the individual pointings last typically for a bit less than 
one day, the above tests will only reveal variability on timescales 
shorter than one day. As a next step we have then performed a 
Eourier analysis of the full (i.e. all four observations) lightcurves 
with 1 ks bins. Eor this purpose we used again the method of 
Heck et al. (119851 1 amended by Cosset et al. (120011 1. To get the 
cleanest results, we again discarded time bins with less than 
80% fractional exposure time. Pig.|9]illustrates the results of this 
Eourier analysis for all four observations. The highest peaks are 
found at frequencies below 0.5 d '. However, the positions of the 
highest peaks and the general shape of the periodograms, espe¬ 
cially above 2d *, are quite different for the three EPIC instru¬ 
ments. As pointed out by Naze et al. (120131 1. having three inde¬ 
pendent cameras onboard XMM-Newton that operate simultane¬ 
ously is extremely beneficial to check the consistency of a timing 
analysis. The differences in Pig.|9] especially at high frequencies, 
most probably reflect different realizations of the noise, whilst 
the low-frequency content of the power spectra hints at the exis¬ 
tence of variations with timescales larger than two days. 

The existence of a long- or medium-term variability can also 
be seen from Table|2] which shows that the mean EPIC count 
rates of A Cep are about 7% lower in observation I than in later 
observations, while a shallower (about 2 - 3%) increase oc¬ 
curs from observation II to IV. The same trends are also seen 
in the fluxes derived from the global spectral fits in Table[3] 
To check whether or not these trends affect the detection of 
shorter periodicities, we have subtracted the average count rate 
from the lightcurve of individual observations and repeated our 
Eourier analysis. The corresponding periodogram is displayed in 
Fig Hni The situation is qualitatively similar to that of Eig.|3 ex¬ 
cept of course for the lack of power at frequencies below about 
0.5 d '. In conclusion, our X-ray dataset alone cannot constrain 
the timescale of the long-term variations and no unambiguous 
short-term periodicity can be identified. 

We can nevertheless fold the X-ray lightcurve with the fre¬ 
quencies detected in our analysis of the simultaneous optical 
spectroscopy. Starting with the vi = 1.315 d * frequency found 
in the variations of our He ii A 4686 data, we note a group of 
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Fig. 9. Fourier analysis of the X-ray lightcurves of A Cep in June 
and August 2013. The top three panels show the power spectra 
for the three EPIC instruments whilst the bottom panel displays 
the power spectral window of the EPIC-pn time series. MOSl 
and MOS2 have very similar spectral windows. The dashed ver¬ 
tical lines indicate the frequencies vi = 1.315, V 2 = 0.245 and 
V 3 = 0.495 d ' found in the analysis of our optical spectra. 






8 10 



(d ‘) 

Fig. 10. Same as Fig.|3 but after subtracting the long-term vari¬ 
ations. 


peaks between 1.1 and 1.4day“^ in Fig.[T0] These peaks are 
probably not significant as their frequencies differ between the 
three EPIC instruments. Pigure[TT]displays the X-ray lightcurve 
for both the observed and the detrendec0 count rates phased with 
vi. The peak-to-peak amplitude of the modulation is 3 - 5 % of 
the EPIC count rate whilst the (one-sided) error bars are 0.8%. 

^ Since the 18.3 hrs period corresponding to vj is shorter than the 
duration of an individual observation, we consider also the results after 
taking out the long-term variations. 


Fig. 11. X-ray lightcurve of the full EPIC (MOSl -H MOS2 H- pn) 
count rate folded in phase with the vi frequency. The data were 
binned into phase intervals of 0.1 and phase 0.0 was arbitrarily 
taken as the starting time of observation I (HJD 2 456 456.1871). 
The bottom and top panels show the results for the background- 
corrected count rates and after removing the long-term trend 
from the background-corrected count rates, respectively (see 
text). 


The amplitude of any modulation of the X-ray flux with vi is 
thus significant at the 2 - 3 cr level only, in agreement with the 
results obtained in our Eourier analysis. 

We have then considered the V 2 = 0.245 d ' frequency which 
we found along with its first overtone in the variations of the Ha 
line profile. Pigure[T2]displays the corresponding lightcurve built 
from the background-corrected count rates. Eor the EPIC count 
rates, the peak-to-peak amplitude now amounts to ~ 9% for a 
mean 0.8% one-sided error bar. This results in a signal signif¬ 
icant at the 5 cr level. This modulation mainly stems from the 
long-term trend discussed above. We have also folded the hard¬ 
ness ratio, computed as the ratio of the EPIC count rates in the 
hard band over those in the soft band, with the V 2 frequency. 
Whilst the source appears harder when the X-ray emission is 
minimium, this result is not significant in view of the large error 
bars. 

Our data sample only 65% of the full 4.1 day cycle, and addi¬ 
tional data are clearly required to establish this periodicity. Still, 
we note that the shape of the EPIC lightcurve in the second panel 
of Pig.[l2] suggests that the modulation could have a double¬ 
wave shape. pjgurefTSl further shows the predicted behaviour of 
the EW of the Ha line. The comparison indicates that the high¬ 
est X-ray flux is observed during the secondary minimum of the 
EW of the optical emission line. This would suggest a close con¬ 
nection, although shifted in phase, with the mechanism respon¬ 
sible for the Ha variability. A similar situation was found by 
Berghofer et al. (119961) for the Ha and soft X-ray emission of 
^ Pup, although more recent X-ray observations did not confirm 
the existence of the 16.7 hrs period. Berghofer et al. (11996b ar¬ 
gued that the phase shift between the two energy domains was 
consistent with the wind flow time between the Ha and X-ray 
emission regions. In the present case, this scenario does not seem 
to work as the X-ray and Ha formation regions overlap to a large 
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0 (^2=0.245 d“‘) 

Fig. 12. X-ray lightcurve of A Cep folded in phase with the V 2 
frequency. The red, blue, green and purple curves correspond 
to the background-corrected count rates as measured during ob¬ 
servations I, II, III and IV respectively. These lightcurves were 
binned into 1 ks bins. The filled symbols with error bars indicate 
the lightcurve binned into phase intervals of 0.1. Phase 0.0 was 
again taken as the starting time of observation I. From bottom to 
top the results are shown for the combination of the two MOS 
detectors, for the pn, for the combination of all three EPIC in¬ 
struments and for the hardness ratio derived for all three EPIC 
cameras. The dashed line in the second panel from the top repre¬ 
sents the variations of the EW of the Ha line scaled to match the 
range of variability of the EPIC count rate and shifted by 0.36 
in phase (with respect to Pig.lHl to match the phase 0.0 of the 
folded X-ray lightcurve. 


extent (see the forthcoming Sect. l5.4l) and wind flow times would 
be much shorter than the observed time lag. Instead, the phase 
shift could rather suggest an anti-correlation between the X-ray 
and Ha emission, e.g. due to an extra heating mechanism. 

In summary, we thus conclude that the X-ray flux of A Cep 
does not undergo highly significant variations on timescales 
of a few hours, but could undergo a rotational modulation on 
timescales of about 4.1 days. This result however calls for con¬ 
firmation through additional X-ray observations. 

A Cep is not the only single O-type star for which a rotational 
modulation of the X-ray emission has possibly been detected. 
Indeed, in addition to the case of ^ Pup (Naze et al. 120131) dis¬ 
cussed above, such a modulation has been claimed in the case 
of ^ Oph (09.5 V) an d of ^ Per (07.5 III(n)((f))). In the former 
case, Oskinova et al. (120011 ) found that the X-ray flux measured 
with ASCA over 1.185 days displays a 20% modulation on a pe¬ 
riod of 0.77 days (i.e. about half the probable rotational period). 
In the latter case, Chandra-HETG data cover a bit less than half 
the probable rotational period of 4.2 days and they display a 15% 
variation. Massa et al. (120141) suggest that the putative rotational 
modulation of ^ Per’s X-ray flux stems from so-called corotat¬ 
ing interaction regions (CIRs). These authors note the absence of 
hardness variations of the X-ray emission of ^ Per and accord¬ 


ingly conclude that the variations of the X-ray flux stem from 
obscuration by optically thick structures. 

5. The X-ray spectrum of A Cep 

In this section, we analyse the X-ray spectrum of A Cep us¬ 
ing techniques of increasing complexity and sophistication. 
Whenever appropriate, we compare our results with those ob¬ 
tained by Bouret et al. (120121) from a combined far-UV, UV and 
optical spectral analysis of A Cep using the CMEGEN model 
atmosphere code (Hillier & Miller [19981) accounting for the ef¬ 
fect of micro-clumping. Bouret et al. (120121) derived an effec¬ 
tive temperature of 36000K, a gravity of log.g = 3.5, a mass- 
loss rate of logM = -5.85 and a wind velocity of 2100kms“'. 
They further determined a moderate chemical enrichment with 
y/yo - 1-33, e(N)/e(N)Q - 8.3 + 3.0. Carbon and oxygen 
were found to be slightly depleted: e(C)/e(C)o = 0.68 + 0.34, 
e(0)/e{0)Q - 0.66 + 0.23, although, within the error bars, these 
numbers are essentially consistent with the corresponding solar 
abundances of Asplund et al. (120051) . 

5.1. The RGS spectrum 

The RGS spectrum (Pig. [13]) provides exploitable data over the 
wavelength range between 9 and 25 A. At the shorter wave¬ 
length end, the flux intrinsically drops off because of the mod¬ 
erate plasma temperatures, leading to too low a S/N ratio in the 
RGS spectrum. Conversely, at wavelengths above 25 A, the X- 
ray emission of A Cep is heavily absorbed by the interstellar 
medium and, to a lesser extent, by the cool material in the stellar 
wind. 

Although rather noisy, the RGS spectrum reveals a num¬ 
ber of emission lines that can be used as diagnostics of the hot 
plasma in the wind of A Cep. Beside several strong Pe xvii lines, 
the strongest features are the Lya lines of N vii, O viii and Ne x, 
and the He-like triplets of O vii, Ne ix and Mg xi. 

The ratio between the fluxes of the hydrogen and helium¬ 
like ions of oxygen and neon can be used to constrain the tem¬ 
perature of the emitting plasma. To do so, we need to correct 
the fluxes for the wavelength dependence of the interstellar and 
wind absorption. Por this purpose, we adopt the interstellar col¬ 
umn density of 0.25 x lO^^cm^^ determined by Bohlin et al. 
( 11978!) and we use the ionized wind column (Naze et al. 12004!) 
of log Awind = 21.79 as derived from the spectral fits in Sect. !5.2! 
The absorption-corrected flux ratios are 3.15^} gg and 0.48^gjg 
for the Ne ix triplet / Ne x Lya and O vii triplet / O viii Lya ra¬ 
tios, respectively. These numbers yield plasma temperatures of 
kT = (0.33 + 0.03) and (0.31 + 0.04) keV for the neon and oxy¬ 
gen ions, thereby indicating that a significant part of the X-ray 
emission of A Cep arises from a plasma at a temperature near 
0.3 keV. 

The strength of the blend between N vii Lya and N vi T 24.9 
compared to the Oviii Lya line further suggests an N/O ratio 
higher than solar. Comparing the emissivities of the nitrogen 
and oxygen lines with the absorption-corrected line strength ra¬ 
tio (0.40;|; j 3 ) as measured on our data, we estimate an N/O ratio 
that would be typically 5.0/[ g times solar if the nitrogen emis¬ 
sion arises from plasma at the same temperature as the oxygen 
lines. This result is somewhat lower than the enrichment found 
by Bouret et al. (!20121 ). but confirms indeed the chemical enrich¬ 
ment of the wind of A Cep. 

The He-like triplets consist of a resonance line (r), an inter¬ 
combination doublet (/) and a forbidden line (/). In the case of 
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Fig. 13. The RGS spectrum of A Cep. The most prominent lines are identified with the labels. The red line shows our best fit with 
the model of Sect. l5.4l The lower panels show the differences between the data and the model. 


A Cep and other massive stars, the / components are strongly 
suppressed, whereas the i components are enhanced compared 
to what one would expect from the transition probabilities. This 
situation stems from the harsh photospheric UV radiation field 
(Gabriel & Jordan 119691 Blumenthal et al. 119721 Porquet et al. 
120011 Leutenegger et al. 12006t that pumps the electrons from the 
upper level of the / transition (2 ^5 1 ) to the upper level of the i 
doublet (2 ^Pi, 2). For instance, the measured "R - f ji ratios of 
the Ovii and Neix triplets are < 0.12 and 0.30 + 0.18, respec¬ 
tively whilst the values expected in the absence of UV pumping 
or collisional depopulation of the 2 \ level would be about 3.8 

and 3 (Porquet et al. 1200111 . 


5.2. Global fits of the X-ray spectrum 

Using the xspec software (Arnaud ri9961 l. we have attempted to 
fit the broadband X-ray spectra of A Cep with apec optically- 
thin thermal plasma models (Smith & Brickhouse 120011) ab¬ 
sorbed by the interstellar medium (with the neutral hydrogen 
column density fixed to 0.25 x 10^^ cm“^) and by an ionized stel¬ 


lar wind using the model of Naze et al. (12004lPl For this purpose, 
we have considered various combinations of the data: 

1. The EPIC spectra of each individual observation. 

2. The EPIC and RGS spectra of each individual observation. 
The RGS spectra were considered over the energy range 0.4 
- 1.5 keV. 

In each case, we tested three different assumptions on the plasma 
abundances: 

1. Solar abundances for all elements following Asplund et al. 
(I2009ll . 

2. Abundances of He, C, N and O as derived by Bouret et al. 
(120121) . all other elements having solar abundances as above. 

3. He abundance fixed to 1.33 times solar following Bouret et 
al. (120121 ) and free N abundance, all other elements having 
solar abundances. 

Prom these tests, we immediately found that the broadband 
spectrum of A Cep requires two plasma components to be rep¬ 
resented by the model. The temperatures of these components 

The stellar wind absorption was imported into xspec as a multi¬ 
plicative tabular model (mtable). 
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Table 3. Spectral parameters for a simultaneous fit of the EPIC and RGS spectra 


Obs. 

log A„ind 

kTi 

norm] 

kTi 

norm2 

{NIH)l(NIHh 

Tv (d-o.f) 

xobs 

Jx 

rcor 

Jx 


(cm-2) 

(keV) 

(cm^*’) 

(keV) 

(cm-°) 



10-" erg cm- 

“ 

I 

21.78);°* 

0.283)°“ 

(2.87);^°) 10-° 

0.86)“ 

(4.12);^°) 10-^ 

3.88)“ 

1.67 (1253) 

6.90 ±0.10 

16.7 

II 

21.81);“ 

0.285);°“ 

(3.68);°°) 10-° 

0.86);°! 

(3.75);!°) lO-'* 

2.91);“ 

1.78 (1265) 

7.24 ± 0.09 

17.8 

III 

21.79);“ 

0.286);°“ 

(3.27)°*) 10-° 

0.86);“ 

(4.07);!°) 10-^ 

3 14'^ 

1.76(1370) 

7.33 ±0.10 

17.8 

IV 

21.79);“ 

u.z,o:7 oo4 

(3.31);°°) 10-° 

0.86);“ 

(4.01);!!) 10-^ 

3 

*^•^^-.60 

1.45 (1279) 

7.45 ±0.10 

18.0 


Notes. The RGS data are fitted between 0.4 and 1.5 keV. The expression of the xspec model is wabs*mtab{wind)*apec(2T). The interstellar 
neutral hydrogen column density is fixed to 0.25 x 10^^ cm“^. The helium abundance by mass y is set to 1.33 solar. Except for nitrogen, all other 
abundances are kept solar. The last two columns provide the observed (last-but-one column) and ISM corrected (last column) fluxes in the energy 

band 0.5 - lOkeV. The normalization of the models is given as 10“''* where and riu are the electron and proton densities of the X-ray 

emitting plasma in cm“^, and d is the distance in cm. The quoted error bars correspond to the 90% confidence level. 
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R„ (R.) 
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Fig. 14. Fits of four individual lines in the RGS spectrum of A Cep. For each row, the top panel illustrates the best-fit model convolved 
with the RGS response and compared with the observed profile. The bottom panel shows the 68.3(1 cr, red continous line), 90 (violet 
dotted line) and 99% (blue dashed line) confidence contours in the t, versus Rq parameter space. The asterisks indicate the best-fit 
parameters. 


are consistently found to he kTi - 0.286 + 0.003 keV for the 
softer component and kT 2 = 0.86 + 0.01 keV for the harder one. 
kT 1 is usually about 0.01 keV lower in the solar abundance case, 
but apart from this small effect, the temperatures are rather inde¬ 
pendent of the assumed composition. They are also independent 
of what observation we consider (see Table[3]l. The best-fit wind 
column density logA^wind = 21.79 + 0.01 is also found to be very 
stable, both with respect to the assumed chemical composition 
and the observation under consideration. Finally, we found that 
the best fit quality was always achieved by the models where 
the He abundance by mass was set to 1.33 times solar and the 
N abundance was left as a free parameter. Its best-fit value was 
found to be 3.3 + 0.4 solar. This is about half the value derived 
by Bouret et al. (I2012I I but is consistent with the value derived 
from the strengths of individual lines (see Sect. EB. Adopting 
the CNO abundances as proposed by these authors resulted in a 
poorer fit of the X-ray spectrum. In this context, we stress that 
the X-ray spectrum of A Cep is not sensitive to the carbon abun¬ 
dance, as the strongest C v and C vi lines are located at wave¬ 
lengths beyond 30 A and fall thus in the heavily absorbed part. 
The same holds for the absorption edges of different ionization 
stages of carbon present in the cool wind. Carbon has therefore 


no measurable impact on the X-ray spectrum of A Cep. However, 
the RGS spectrum is quite sensitive to oxygen. If we let the abun¬ 
dance of oxygen vary in the fit, we obtain an abundance slightly 
higher than solar, but fully consistent with the solar abundance 
within the errors. In Table[3] we present results with an oxygen 
abundance set to solar. 

From the absorption-corrected X-ray fluxes in Table^] and 
a distance of 0.95 kpc (as adopted by Bouret et al. I2012rh . one 
can estimate an X-ray luminosity of 1.9 x 10^^ erg s“' in the 0.5 
- lOkeV energy band. Bouret et al. (120121 1 infer a bolometric 
luminosity of log ^ = 5.80. This in turn leads to log = 
-7.10 + 0.02 for A Cep. Note that the quoted error only accounts 
for the dispersion in the value of the absorption-corrected X-ray 
fluxes. It does not account for any uncertainties on the bolomet¬ 
ric luminosity or the stellar distance. This ratio is towards the 
fainter end of, but still compatible with, the average relation for 
O-stars observed with XMM-Newton as reported by Naze (120091 1. 
Therefore, the level of X-ray emission of A Cep appears rather 
normal for a single O-star, and there is no overluminosity due 


^ This is significantly larger than the Hipparcos distance of (606 ± 
81) pc (van Leeuwen l2007t . 
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to a magnetically confined wind, contrary to what is seen in a 
sample of OB stars with magnetic detections (Naze et al. 120141 1. 

5.3. Individual line profiles 

As a next step in the analysis of the RGS spectrum of A Cep, 
we have performed fits of individual spectral lines. If the X-ray 
emission arises from a hot plasma distributed throughout the oth¬ 
erwise cool wind and moving along with the latter, we can ex¬ 
pect the spectral lines to be broad and skewed (MacFarlane et 
al. 119911 Owocki & Cohen 12001112006l l. Fitting the lines with 
such a model can provide us with a first order estimate of where 
the X-ray plasma is locatec0, and most of all, it can show us 
whether or not the lines are indeed broad and skewed. Given 
the limited quality of our data, we have fitted the profiles with 
the simplest possible model. Previous studies (Sundqvist et al. 
120121 Herve et al. 120121 1 have shown that, unless the porosity 
lengths are very large, porosity has a rather marginal effect on X- 
ray line profiles. Large porosity lengths are ruled out by the ob¬ 
served line profiles (Herve et al. 120131 Leutenegger et al. l2013l l. 
In our calculations, we thus assume that there is no porosity. We 
further assume that the plasma filling factor has no radial depen¬ 
dence. The velocity of both the hot plasma and the cool wind 
is taken as v{r) - Vco(l - with Voo = 2100kms“' (Bouret 
et al. 120121 1. The only remaining free parameters are the inner 
radius of the emission region Rq and the typical optical depth 
of the wind t, = We have fitted the profiles using our 

own code where the synthetic profiles are convolved with the 
line-spread function of the RGS (FWHM = 0.07 A) and are then 
compared to the observed profile. As a consistency check, we 
have also used the windprofile modelfl implemented under 
xspec. Both methods give generally consistent results. The re¬ 
sults of some of our fits are illustrated in Figs.fT4l-fT6l 

One of the strongest lines in the RGS spectrum of A Cep is 
the Fexvii blend near 15 A. Whilst this complex is dominated 
by the blend of the Fexvii AA 15.014, 15.265 lines, other tran¬ 
sitions, such as Oviii Lyy at 15.175 A probably also contribute. 
This makes this feature not well suited for line profile fitting. 
The same restriction applies to the Ne ix He-like triplet which is 
blended with Fe xvii lines. 

The Fexvii lines around 17 A are in principle better candi¬ 
dates for line profile fitting. The Fexvii A 16.78 line is relatively 
well isolated. On the other hand, Fexvii AA 17.051, 17.096 are 
heavily blended. The latter of these lines arises from a metastable 
level that can be depleted in the presence of a strong photo- 
spheric UV radiation field (Mauche et al. 120011 Herve et al. 
I2013I I. As a result, the intensity of the A 17.096 line relative to 
T 17.051 is difficult to predict, making this blend more prob¬ 
lematic for line profile fitting. Comparing the ratio between 
the observed fluxes of the Fexvii AA 17.051, 17.096 blend and 
the Fexvii d 16.78 line to the corresponding ratio of intrinsic 
emissivities of these features, we empirically estimate that the 
A 17.096 line has its contribution reduced by at least a factor two 
with respect to its theoretical emissivity. Fig. [14] illustrates the 
fits obtained under this assumption. 

The three components of the O vii He-like triplet were fitted 
simultaneously assuming identical line parameters (t,, Rq) for 
all three lines. As can be seen on Fig.[T5| the forbidden line is 


® Although we caution that only a fit of the full spectrum consistently 
deals with the blends and contributions from plasmas at different tem¬ 
peratures and positions in the wind, see Herve et al. ( 120131 and Sect. l5.4l 
’ http://heasarc.gsfc.nasa.gov/xanadu/xspec/models/windprof.html 
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Fig. 15. Same as Fig.[T4| but for the He-like triplet of O vii. 


strongly suppressed with a Icr upper limit on the R — fli ratio of 
0.12. Accounting for the UV pumping by a diluted photospheric 
radiation field, the upper limit on R. translates into an upper limit 
on Rq of 12.5 Rt, which is entirely consistent with the best-fit 
Rq - 2.85 R, derived from the fit of the line profile. 

The situation is more complicated in the case of the N vii Lya 
line which is blended with N vi 424.9. If we assume that both 
lines arise in the same region of the wind, we can estimate the 
contribution of the N vi 4 24.9 line to the blend from the temper¬ 
ature dependence of the emissivities derived in our global fit of 
the broadband spectra (Sect. l5.2l l or in the global fit of the RGS 
spectrum (see the forthcoming Sect. l5.4b . In the former case, the 
combined EPIC and RGS spectra were fitted with plasma com¬ 
ponents at temperatures of 0.29 and 0.86keV. At these tempera¬ 
tures, the contribution of the N vi line is expected to be respec¬ 
tively ~ 60 and 400 times lower than that of the Lycr line. If we 
consider instead the plasma temperatures found in our global fit 
of the RGS spectrum (Sect. l5.4T l. we expect the strongest con¬ 
tribution of the N VI line to arise from the 0.20 keV component. 
However, even in this case, the emissivity of the N vii Lya dou¬ 
blet exceeds that of the N vi 424.9 by a factor ~ 20. In the case 
of 4 Cep, the N vi 4 24.9 line should thus play a negligible role, 
suggesting that we can treat this blend as consisting only of the 
N VII Lya doublet. The results obtained under this assumption 
are shown in Fig. [16] As one can see on this figure, the fit is of 
quite limited quality, especially in the blue part. The best-fit pa¬ 
rameters differ significantly from those of the other lines. The 
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Fig. 16. Best fit for the N vii Lya line, assuming no contamina¬ 
tion by N VI /124.9. 



Fig. 17. Parameters of the best fits for those lines that could be 
fitted. The continuous line in the lower panel illustrates the typ¬ 
ical optical depth as computed from the CMFGEN model used 
by Bouret et al. (120121 1 and adopting an opacity k evaluated at a 
distance of 2 R* from the centre of the star. 


best-fit wind optical depth is now 0.0 and the best-fit inner ra¬ 
dius of the emitting region is no longer constrained, except that 
it must be larger than 8.5 R*. 

To improve the fit of the blue part of the profile and simulta¬ 
neously bring the fit parameters into better agreement with the 
trends seen in the other lines, we would need to assume a sig¬ 
nificantly larger contribution of the N vi A2A.9 line. As pointed 
out above, the known temperature components do not allow 
for such a contribution. Yet, it could arise from a significantly 
cooler plasma component that would otherwise be hidden by the 
large interstellar absorption that prevents us from seeing the N vi 
He-like triplet emission. In the case of ^ Pup, Herve et al. (I2013I I 
found that a fit of the RGS spectrum required a four-temperature 
plasma, with the coolest component at O.lOkeV accounting for 
the bulk of the N vi He-like triplet emission. At this plasma 
temperature, the contribution of N vi T 24.9 would be about half 
that of the N vii Lya doublet. To investigate the possibility that 
such a cool component might also exist in the wind of A Cep, 
we have repeated the fits of the N vii Lya -i- N vi T 24.9 blend 
assuming a relative contribution of the N vi feature varying from 
0.0 to 0.9 times the contribution of N vii. Increasing the relative 


contribution of N vi leads to an increase of t* as expected, but 
Rq remains at the upper limit of our parameter grid and the fit 
quality actually degrades. Hence, we conclude that the present 
model cannot easily explain the morphology of this line un¬ 
less its parameters deviate significantly from those of other lines. 


Table 4. Best-fit parameters for the lines where a good fit quality 
was achieved. 


Line Rq (R.) t. 


Nex 4 12.134 

1 05 

^•^^-0.00 

2.60),,0 

O VIII A 16.006 

1 25‘^" 
^•^^-0.20 


Fexvii A 16.780 

^•^'^-1.00 

^•^^-1.25 

Fexvii 44 17.051-17.096 

^•^'^-1.00 

^•^'^-0.85 

Oviii 4 18.970 

^•^*^-0.95 

2 -20^uo 

O VII 44 21.602-22.098 

9 

1.40-,3 

N VII 4 24.780 

8.85+- 

0.00 o.oo 


Notes. Rq is measured from the centre of the star. Whenever possible, 
the Icr errors are estimated from the contours in the parameter space. A 
‘...’ indicates situations where the Icr contours extend beyond the range 
of parameters explored in our calculations. 


TableH] summarizes the results of our fits and Lig.[T7] illus¬ 
trates the overall behaviour of the best-fit t, and Rq values as 
a function of wavelength for those lines where a good fit was 
achieved. Most transitions have best-fit Rq values consistent with 
an emission region starting between 2 and 3 R*. The only excep¬ 
tions are the shortest wavelength transitions that seem to orig¬ 
inate very near the photosphere, although, also in these cases, 
the 1 cr confidence ranges are still consistent with an onset radius 
between 2 and 3 R, from the centre of the star. There is no obvi¬ 
ous trend of t» with wavelength, though we note that the (large) 
error bars overlap with the expected t,(/1 ) relation at a distance 
of 2 R* from the centre of the star assuming the mass-loss rate 
inferred by Bouret et al. (120121) . 

5 . 4 . Global fit of the RGS spectrum 

Under the assumption of a wind embedded plasma, the X-ray 
spectrum of a single non-magnetic O-type star stems from the 
combination of emissions arising from plasma at different tem¬ 
peratures and different locations in the wind. Hence, the above 
approach of fitting individual line profiles only provides a first 
order approximation of the plasma properties. A more rigorous 
approach is therefore needed. Herve et al. (12012112013112015t 
have developed a code that aims at achieving a consistent fit of 
the high-resolution X-ray spectrum, simultaneously accounting 
for the presence of several plasma components with different lo¬ 
cations inside the wind. In this section, we apply this code to the 
specific case of the RGS spectrum of A Cep. 

In our model, the part of the wind that contains the X- 
ray plasma is divided into a number of spherically symmetric 
shells. To allow for a position-dependent filling factoi0 of the 
hot plasma, we adopt shells with a minimum width of 0.1 R, 
(see also below). The X-ray emissivities of each of these shells 
are computed assuming an optically thin thermal plasma (using 
the AtomDB code, see Loster et al. l2012l and references therein). 


* The filling factor is defined here as the ratio of the X-ray emitting 
volume and the total volume of the shell under consideration. 
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The important parameters are the plasma temperature and the 
quantity of hot gas contained in the shell. The shells are assumed 
to move along with the ambient cool wind and their emission 
is Doppler-shifted into the frame of reference of the observer. 
The optical depth of the cool wind material is computed along 
each line of sight from the observer to an emitting shell using 
the k(A, r) relations computed with CMFGEN (see Herve et al. 
l20T2ll20T^l20T5l for details). 

The UV pumping of the 2 i level in He-like ions to their 
2 level and the resulting enhancement of the intercombina¬ 
tion line at the expense of the forbidden line are taken into ac¬ 
count following the approach of Leutenegger et al. (120061 1. For 
this purpose, we adopt the UV fluxes from our CMFGEN model. 

The RGS spectrum features lines that arise from plasmas at 
different temperatures (e.g. Fe xvii lines are preferentially pro¬ 
duced by a hot plasma with kT close to 0.6 - 0.7 keV, whereas 
the emissivity of N vii lines peaks around 0.2 keV). To reproduce 
the observed spectrum, we therefore need to combine models 
with different temperatures and with different radial dependen¬ 
cies of their Ailing factors. We thus need to find a combination of 
our models that produces a synthetic spectrum that best matches 
our observed spectrum. To this aim, we have to build a database 
which samples the different parameters of the synthetic spectra 
simulator code such as the inner and outer boundary of the emis¬ 
sion region, the abundances of CNO, the mass-loss rate, the wind 
velocity, the porosity length and, of course, the X-ray emitting 
plasma temperature. In the present case, several parameters were 
frozen. Indeed, we adopted the wind velocity law of Bouret et 
al. (120121) and we neglected porosity (see Sect. 15.31 ). Solar CNO 
abundances did not allow to fit the data, and we thus adopted the 
same chemical composition as derived by Bouret et al. (12012t . 
Our fitting routine, based on the Bound Variable Least Square al¬ 
gorithm, systematically explores this database, combining up to 
four plasma components with the same values for the chemical 
composition and mass-loss rate. 



r (R.) 


Fig. 18. Radial distribution of the Ailing factor of the hot plasma 
inside the stellar wind of A Cep as inferi'ed from our fit of the 
RGS spectrum shown in Fig.fTSl 

Overall our best-fit model does a good job in reproducing 
the observed RGS spectrum (see Fig.fTJt. Still it fails to re- 



1B.8 18.9 19 19.1 


X (A) 

Fig. 19. Impact of some model assumptions on the fits of the 
RGS spectrum of A Cep, illustrated by the example of the O viii 
Lya line. In the top panel, the red solid line illustrates the best- 
fit model (assuming minimum shell width of 0.1 R,), whilst the 
blue dashed line corresponds to the best fit for a minimum shell 
width of IR,. The histogram yields the observed line profile. 
The bottom panel shows the same figure, except for the blue 
dashed line which corresponds now to the best fit where we have 
artificially set the onset radius of the X-ray emission region to 
be located at r > 2 R,. 


produce some features. For instance, whilst the strong Fexvii 
/115.014 line is very well reproduced, the strengths of the Fe xvii 
AA 15.265, 16.780, 17.051 lines are overpredicted by our model. 
For the Fexvii A 17.051 line, we found that the contribution of 
Fe XVII A 17.096 must be essentially negligible as a result of UV 
pumping of the upper level of this transition. 

The best-fit mass-loss rate of A Cep was found to be 1.5 x 
10 “® Mq yr“', in excellent agreement with previous studies (1.5x 
10 -®Moyr-', Sundqvist et al. l20Tl1lm^ 1.4 x lO^^Moyr-^, 
Bouret et al. 120121 1.6 x lO^^Moyr ', Surlan et al. 120131 ) that 
deal with the dumpiness of the stellar wind in different wayfl 
This agreement indicates that high-resolution X-ray spectra can 
indeed be used to constrain the mass-loss rates of O-type stars. 
The above values of M are lower than the theoretical mass-loss 
rate based on the recipe of Vink et al. (120001 ). Indeed, applying 
the latter recipe to the stellar parameters as inferred by Bouret 
et al. (12012t yields 4.5 x 10“^ Mq yr“', i.e. a factor 3 larger than 
what we have found here. Various studies comparing mass-loss 
rates derived from observations with those predicted by the Vink 
et al. (120001) formalism have revealed similar discrepancies (e.g. 
Bouret et al. 120121 Puls et al. 120151 ). with the Vink et al. (120001 ) 
rates being a factor 2-3 too large. Our results indicate that A Cep 
is yet another case where this conclusion applies. 

^ Sundqvist et al. ( 12011112014b and Surlan et al. ( 120131 indepen¬ 
dently analysed UV wind resonance lines accounting for the effects 
of clumps of any optical depth (from optically thin to optically thick). 
Sundqvist et al. ( l2011ll20T4l also account for porosity in velocity space, 
whilst Surlan et al. ( 120131 considered the velocity dispersion inside the 
clumps. Bouret et al. ( 120121 simultaneously fitted the FUV, UV and op¬ 
tical spectra adopting a parametric filling-factor approach and assuming 
optically thin clumps. 
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Three plasma components were needed to achieve the fit 
shown in Fig.[T3] The best-fit plasma temperatures were 0.200, 
0.272 and 0.684keV. The radial distributions of their hlling fac¬ 
tors are shown in Fig. [18] Although we caution that this distri¬ 
bution is likely not unique (see below), it is worth pointing out 
that the bulk of the hot plasma is located within less than 2 R, 
above the stellar surface and some emission can actually arise 
very close to the stellar surface. The same conclusion is reached 
for models which have slightly different temperatures but pro¬ 
duce fits of almost equal quality. This result also agrees quite 
well with the Rq values found in Fig. [n] 

Although our algorithm systematically explores the parame¬ 
ter space to search for the best combination of a predefined num¬ 
ber of plasma components, we caution that the distribution in 
Fig. (HI is likely not unique. For instance, while we found that 
two plasma temperatures are clearly not sufficient to reproduce 
the observed spectrum, the possible presence of a fourth lower- 
temperature component cannot be constrained. Moreover, the 
limited S/N ratio of the RGS spectrum of A Cep and the strong 
interstellar absorption prevent us from gathering stronger con¬ 
straints on the temperatures and locations of the emitting plasma 
that would lift some degeneracies on the radial distribution of 
the filling factor. The latter might also depend on some of our 
assumptions, such as the minimum width of the emitting shells. 
Indeed, Fig. [18] shows that the various emission regions extend 
only over a few times the minimum shell width. To quantify 
the sensitivity of our fit on the assumed minimum width of the 
shells, we have performed another ht setting the minimum shell 
width artificially to 1 R,. The top panel of Fig. [19] illustrates the 
comparison between the best fit obtained under this assumption 
and the one of Fig. [13] To enhance the visibility of this compar¬ 
ison, we focus on a single line (O viii Lya), although we stress 
that the full spectrum was fitted. As can be seen, except for a 
marginal change in line width, there is essentially no difference 
between the line profiles obtained for both assumptions. This is 
actually not surprising as Fig.[T8]reveals shells of the dominant 
plasma component (the one at 0.272 keV) to be present between 
r = 1.1 and r = l.OR*, i.e. over a range of almost 1 R*. We thus 
conclude that the detailed radial distribution of the hot plasma 
hlling factor depends on the assumptions that we make on the 
shell width. 

How robust is then the hnding that the X-ray emission starts 
already at 1.1 R,? To check this, we have done yet another ht 
where we artihcially require the onset radius of the X-ray emis¬ 
sion region to be located at r > 2 R*. The result is shown in the 
bottom panel of Fig. [19] This time, we see that a larger onset ra¬ 
dius leads to a blue-shift of the line centroid. Whilst this effect is 
small if we consider a single line (as shown in Fig.[T9]). consid¬ 
ering the spectrum as a whole leads to a better ht for the onset 
radius at 1.1 R*. 

At hrst sight, the onset of the X-ray emission that close to 
the stellar surface is a challenge for theoretical models. Whilst 
the hydrodynamical simulations of Feldmeier et al. (119971) pre¬ 
dict that the instability of the line driving mechanism already 
produces considerable structures as close as 0.1 R, above the 
photosphere, in these models, the bulk of the X-ray emission 
arises between 1 and 3 R, above the photosphere. More re¬ 
cently, Sundqvist & Owocki (120131 ) included the effect of limb- 
darkening in their hydrodynamical simulations. These authors 
noted that limb darkening favours the development of wind 
structures already at r = 1.1 R*, although these simulations pro¬ 
vide no information about the level of X-ray emission. Such an 
early onset of the X-ray emission could also affect the formation 
of important UV resonance lines as well as optical lines formed 


at the base of the wind. From the observational point of view, 
we note that similar situations have been reported for other stars 
using independent methods. For instance, in their analysis of the 
09.5 II star 6 Ori, Shenar et al. (120151) derived an onset radius of 
~ I.IR.. 

Finally, one may wonder about the signihcance of the outer 
boundary of the X-ray emitting plasma near 2.5 R,. In some 
models, we found indeed plasma components that are present 
out to much larger radii. However, in these cases, the corre¬ 
sponding hlling factors were very low so that these regions con¬ 
tributed only a marginal fraction of the line Hux. Leutenegger 
et al. (120131 ) showed that the existence of an upper radial cutoff 
impacts the shape of the lines (their Fig. 7), leading to a more 
asymmetric prohle for a given value of t*. For our analysis of 
Sect. 15.31 where we made the assumption of an inhnite emission 
region, this may imply an overestimate of the values of t,. 

In summary, we thus conclude that the results that the X- 
ray emission of A Cep starts already at 1.1 R, and arises from a 
region that spans about 2 R, in the wind are quite robust. 

This distribution is actually quite different from the one ob¬ 
tained in the case of ^ Pup (Herve et al. 12013t . In the latter star, 
the softer emission starts at 0.5 R« above the photosphere and 
reaches out to very large distances in the wind. The harder emis¬ 
sion component in the spectrum of ^ Pup arises from between 
2 and 3 R* above the photosphere. We thus see that the X-ray 
emitting plasma is much closer to the stellar surface in the case 
of A Cep and covers a narrower range of distances than in the 
wind of ^ Pup. 

6. Conclusion 

In this paper, we have presented the hrst-ever X-ray analysis of 
the 06 lef star A Cep. Because of its spectral type, this star is 
certainly the closest cousin of ^ Pup, although, because of its 
larger distance and higher interstellar absorption, its X-ray flux 
at Earth is much lower. We have further gathered coordinated op¬ 
tical spectroscopy to monitor the variability of A Cep during our 
fom XMM-Newton observations. These data led to the following 
results: 

- Optical spectroscopy revealed variability of the Ha and He ii 
/14686 emission lines. Whilst the former appears to vary on 
a timescale of 4.1 days, very likely associated with the ro¬ 
tation of the star, the latter varies on an unrelated timescale 
(18.3 hrs). Given the proximity of the formation regions of 
both lines, as derived from 1-D model atmosphere codes, 
this difference is a puzzle that needs to be addressed in future 
studies. 

- While there is only marginal X-ray variability during in¬ 
dividual X-ray pointings, we find a signihcant (5cr) inter¬ 
pointing variability. Although the timescale of these varia¬ 
tions cannot be constrained by the sole X-ray data, it turns 
out to be compatible with the 4.1 days timescale derived 
from the Ha data. One of the maxima of the X-ray emission 
seems to occur when the EW of Ha is minimum, indicating 
a phase shift between the two variations. 

- The high-resolution X-ray spectrum of A Cep reveals broad 
and slightly skewed emission lines. As in the case of ( Pup, 
most of these lines are rather well reproduced by a model that 
assumes to first order a spherically symmetric distribution of 
shock-heated plasma embedded in the cold wind and moving 
along with the latter. 

- A global ht of the RGS spectrum indicates that the bulk of 
the hot plasma must be located close to the stellar surface and 
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yields a mass-loss rate in very good agreement with previous 

determinations based on optical and UV spectra. 

At first sight, a rotational modulation of the X-ray flux could 
hint at the effect of wind-conhnement by a moderately strong 
dipolar magnetic field (Babel & Montmerle 119971 ud-Doula & 
Owocki l2002t which would force the wind, out to the Alfven ra¬ 
dius, into rigid co-rotation with the star. However, such a global 
magnetic held has not been detected (e.g. David-Uraz et al. 
12014b . Furthermore, if this scenario applies in the case of Oef 
stars, the magnetically conhned wind does not extend over a 
wide part of the X-ray emission region. Otherwise, one would 
expect to observe a strong X-ray overluminosity (Naze et al. 
12014b along with some narrow, symmetric lines, as seen notably 
for 0^ OriC (Gagne et al. [2005ll and HD 148 937 (Naze et al. 
I2012bb . However, this is dehnitely not the case neither for C, Pup 
(Naze et al. l2012al Herve et al. 12013b . nor for A Cep, as we have 
shown here. 



Fig. 20. Cartoon of the CIRs in the wind of A Cep and its impact 
on the X-ray and Ha emission. The Ha emission region is shown 
by the blue hatched area. Outside the CIR, the X-ray emission 
region overlaps with the Ha emission region. Near the base of 
the CIR, excess X-ray emission is produced and Ha emission is 
reduced. 

Henrichs & Sudnik (12014b argue that localized magnetic 
helds might exist in the case of A Cep, leading to spots corotat¬ 
ing with the star that might be at the origin of corotating interac¬ 
tion regions (Cranmer & Owocki fT996l Lobel & Blomme l2008b . 
Harries ( 12000b performed Monte-Carlo line prohle calculations 
for a single corotating spiral density enhancement in the wind of 
an O supergiant. He showed that such a one-armed density wave 
actually produces a double-wave modulation of the equivalent 
width similar to what we have observed for the Ha line of A Cep. 
In the simulated line prohles, minimum EW occurs at phases at 
which the base of the spiral is eclipsed by the star or crosses in 
front of the latter (Harries r2000b . Of course, such double-wave 
modulations could also arise if there were two spots. Indeed, 
Lobel & Blomme (12008b performed 3-D radiative transfer calcu¬ 
lations with hydrodynamical models to simulate the discrete ab¬ 
sorption components due to CIRs in the BO.5 Ib star HD 64760. 


These authors tested models with one or two spots (located on 
opposite sides of the stellar equator) and found a slightly better 
agreement between their simulations and the observed variations 
for the two-spot modeQ 

The question that arises is whether the same CIR could im¬ 
pact the level of the observed X-ray emission. Our results sug¬ 
gest that this is the case. Such a modulation could be either a 
consequence of extra emission produced in the CIR itself or of 
the associated density enhancement which would lead to a mod¬ 
ulation of the column-density towards the hot plasma. A can¬ 
didate for the first scenario could be X-ray emission produced 
behind the shock associated with the CIR. Heating the plasma 
to kT - 0.200 via this mechanism requires a velocity jump of 
400 km s ^ Such values were deemed possible in the work of 
Cranmer & Owocki (11996b . However, at least for HD 64 760, the 
more recent hydrodynamical simulations of Lobel & Blomme 
(12008b yield maximum velocity differences across the CIR that 
do not exceed 140km s '. Such values would not be sufhcient to 
provide a signihcant X-ray emission. An alternative scenario to 
produce the additional X-ray emission could be the action of the 
localized magnetic helds that were proposed to produce the spots 
that generate the CIRs (Henrichs & Sudnik 120 14b . An excess of 
magnetically-heated plasma near the base of the CIR could al¬ 
ter the local ionization, thereby creating a hole in the Ha emis¬ 
sion region which could explain the apparent anti-correlation be¬ 
tween the X-ray and Ha emissions. 

Ligurel20l illustrates a possible conhguration. The shape of 
the CIR in Lig.|20]is based on the expression of the spiral streak¬ 
line of Lullerton et al. (11997b . In Sect. l5.4l we found that the bulk 
of the X-ray emission of A Cep arises from the inner parts of the 
wind. Therefore, in view of Lig.|20l it appears unlikely that the 
X-ray modulation could be due to absorption by the density en¬ 
hancement associated with the curved part of the CIR. Indeed, 
maximum X-ray emission is observed when the base of the CIR 
is roughly aligned with the line of sight and the projection of 
the curved part of the CIR on the plane perpendicular to the line 
of sight covers a wide area. A more plausible explanation would 
thus be that the base of the CIR produces itself some X-ray emis¬ 
sion. 

Considering the results of Naze et al. (12013b on ^ Pup, Massa 
et al. (12014b on ^ Per, and our own results on A Cep, obser¬ 
vational evidence is accumulating for the impact of large-scale 
wind structures on the X-ray emission of O-type stars even 
though current theoretical models have difficulties to explain 
this. Because such structures might be transient, coordinated 
optical and X-ray observations are mandatory to hrmly estab¬ 
lish this link. Such studies are pushing the XMM-Newton and 
Chandra observatories to the limits of their capabilities, but are 
clearly worth the effort. In the more distant future, these studies 
will greatly beneht from the unprecedented collecting area of 
ESA’s Athena X-ray observatory which will allow us to observe 
the impact of the structures on individual X-ray emission lines 
(see Sciortino et al. 12013b . 
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